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Incident  wave  amplitude 
Gravitational  acceleration 
Body  draft 

Body  moment  of  inertia  in  pitch  about  the  center  of  gravity 

=  n/TT  \  „ 

Bessel  function  of  the  first  kind  of  order  zero 

’  •  2  • 

Wave" number,  w  /g 

2 

Radius  of  gyration,  =  I/m 
Body  mass 

Unit  normal  vector  into  the  body 

(z  -  S(z)  dz 

-H  ^ 
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p  Pressure 

Qj^(k)  =  (z  -  Zq)^  S(z)e^^  dz 


R(z)  Sectional  radius  of  the  body 

S(z)  Sectional  area  of  the  body 

r  Polar  radius,  4-  y^ 

t  Time 

(x,  y,z)  Cartesian  coordinate  system 

Zq  Coordinate  of  the  center  of  gravity 

^  Heave  displacement 

C*  free  surface  elevation 

0  Polar  coordinate 

t  Surge  displacement 

p  Fluid  density 

‘I’  Velocity  potential 

X  Vertical  prismatic  coefficient 

([)  Pitch  angle 

(jj  I'T'oquency  of  oscillations 
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A  linearized  theory  is  developed  for  the  motions  of  a 
slender  body  of  revolution,  with  vertical  axis,  which  is  float¬ 
ing  in  the  presence  of  regular  waves.  Equations  of  motion 
are  derived  which  are  undamped  to  first  order  in  the  body 
diameter,  but  second-order  damping  forces  are  derived  to 
provide  solutions  valid  at  all  frequencies  including  resonance. 
Calculations  made  for  a  particular  circular  cylinder  show 
extremely  stable  motions  except  for  the  low  frequency  range 
where  very  sharp  maxima  occur  at  resonance. 

o 

INTRODUCTION  ' 

The  motions  of  a  vertical  body  of  revolution,  which  is  floating  in  the 
presence  of  waves,  present  a  problem  of  interest  in  several  connections". 
The  motions  of  a  spar  buoy,  of  a  wave-height  polo,  and  of  floating  rocket 
vehicles  are  important  examples  of  such  a  problem.  The  same  methods 
developed  for  these  motions  may  be  applied  to  find  the  forces  acting  on 
offshore  radar  and  oil-dri.’ling  structures. 

A  theoretical  di.scu8sion  of  this  problem,  which  also  treats  the  sta¬ 
tistical  problem  of  motions  in  irregular  waves,  has  been  presented  by 
Barakat.*  However,  this  ana.lysis  is  rcostricted  to'the  case  of  a  circular 

0  O 

cylinder  find  is  based  upon  .several  semi -empirical  concepts  of  applied 
ship-motion  theory.  An  .alternative  proc  edure  is  to  formulate  the  (inviscid) 
hydrodynamic  problem  as  a  boundary- value  problem  for  the  velocity  po¬ 
tential  and  to  employ  slender-body  techniques  to  solve  this  problem.  The 
latter  approach  is  followed  in  the  present  work,  leading  to  linearized  equa¬ 
tions  of  motion  which  may  be  solved  for  an  arbitrary  slender  body  with  a 
vertical  axis  of  rotational  symmetry.  The  particular  case  of  a  circular 
cylinder,  whose  centers  of  buoyancy  and  gravity  coincide,  is  treated  in 
detail  and  curvc.s  are  presented  for  the  amplitudes  of  surge,  heave,  and 
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pitch  oscillations. 

In  deriving  the  hydrodynamic  forces  and  moments  acting  on  the  body, 
we  shall  assume  that  the  incident  waves  and  the  oscillations  of  the  body  are 
small,  and  thus  we  shall  retain  only  terms  of  first  order  in  these  ampli¬ 
tudes.  We  shall  also  assume  that  the  body  is  jslender.  The  analysis  with 
only  first-order  terms  in  the  body^s  diameter  leads  to  undamped  resonance 
oscillations  of  infinite  amplitude.  To  analyze  the  motions  near  resonance, 
it  is  necessary  to  introduce  damping  forges  which  are  of  second  order  with 
respect  to  the  diameter  "length  ratio. 


THE,  FlRST"ORDER  VEEQCITY 

We  shall  consider  the  hydrody^namie,  problem  of  a  floating  slender 
body  of  revolution  with  a  'V’erticai  aKss  in  the  presence  of  small  incident 
surface  waves..  Let  (x^y^zj  be  a  fixed  Cartesian  Coordinate  system  with 
the  z-axis  positive  upwards  and  the  plane  z  u  0  situated  at  the  undisturbed 
level  of  the  free  surface.  The  X'>axis  is  taken  to  be  the  direction  of  propa'-> 
gation  of  the  incident  wave  system,  and  motion  ®f  the  body  is  assurned 
t®  be  confined  tv  the  plane  y  “  0.  We  shall  also  em,ploy  a  co©rdii»ate  sys« 
«em  fx^y^f  fixed  in  the  body,  with  z”  the  axis  of  the  body,  s®  that  with 
the  body  at  rest,  fx,y,z)  =  fx'.yhz');  and  a  circular  cylindrical  system 
Ci'sO.z},  where  x  =  r  cos  0  and  y  ~  r  sin  0,  If  £,  £;,  and  are  the  instant 
'tane«n»s  aanplitudes  «>f  surge,  heave,  and  pitch,  respectively,,  rclatsveoto' 
the  body's  center  of  gravity,  it  follows  that 

"  x;  JP,  -f-  x'  <ttss  njf  -S'  zq)  sin  ifs 

ofay*  -HI 

O  O 

a  2  ^  e,  54®  sisb  4^  +  -  Z%3  COS  4*  «  ^  Oo 

o  °  “ 

wheS’e  Zq  ss^the  vertica!  conrdtda^e  (pf  the  Center  of  gravity  in  the  bodys 
fvxed  system;  see  F^igUre  5„  The  (displacements  and  i|j  are  assumed 

t®  be  small  oscillatory  fanstiwrts  cf  time;  Wc  shall  consistently  linearize 
by  neglegting  terms  of  second  order  in  these  functions  or  their  products 
with  the  incident  wave  amplitude  Aa  Thus  Sqwatioit  1 1|  onay  be  replaced  by 

o  ^  n 

o 

’’  ^ 

o 

T  O 

O 


0  fi 


3c^  |  +  +  Zq) 

r^Y"  [2] 

z  =  ^  —  xP  +  z' 

o 

I:?  an^deal  ingomp^Fessoble  fiujd  is  assurjied.,  there  exists  a  velocity 

C) 

potential,  $ (x<oy,  z,  tj,,  satisfying  Laplade^s  equation,  su(?h  that  its  gradi¬ 
ent  is  equal  ta  the  velocity  of  the  fl’taid,,  This  function  must  satisfy  the 
following  boundary  Conditions^  .) 

fll  On  the  body,  the  normal  velocity*  component  of  the  body  must  equal 
the  normal  derivative  of  For  a  body  of  revolution  defined  by  the  equa¬ 
tion  r'  -R(z*),  where  r®  -  Vx’'^  +  y*'',  this  boundary  condition  may  be 
expressed  by  the  equation*  ” 
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[3] 


D 

Dt 


R(z')J 


A  +  V)[r'  -  R(z’)]  -  0 

on  r'  =  R(z') 


(2)  On  the  free  surface,  the  normal  velocity  component  of  the  free 
surface  must  equal  the  normal  velocity  component  of  the  fluid  particles 
in  this  surface,  and  the  pressure  must  equal  atmospheric  pressure.  In 
the  linearized  theory,  these  conditions  reduce^  to 


-  +  G  - 

at2 


0  on  z  =  0 , 


[4] 


or  in  the  case  of  a  sinusoidal  disturbance  with  frequency  u>, 


K  $  -  — *  =  0onz=0,  [5] 

oz 

where  K  =  oj  /g. 

(3)  At  infinite  distance  from  the  body,  the  waves  generated  by  the  body 
are  outgoing  (the  radiation  condition). 

The  free  surface  condition,  Equation  [5],  and  the  radiation  condition 
are  satisfied  by  the  potential  of  oscillating  singularities  beneath  the  free 
surface;  the  boundary  condition  on  the  body  may  be  satisfied  by  a  proper 
distribution  of  these  singularities.  This  distribution  may  be  found  from 
slender-body  theory  but  some  care  is  required  in  linearizing  the  present 
problem.  If  r'  =  R{z')  is  the  equation  of  the  body  surface  over  its  sub¬ 
merged  length  (-H  <  z'  <  0),  we  shall  assume  that  R  and  its  first  deriva¬ 
tive  are  continuous,  that  R{  — H)  =  0,  and  that  the  magnitude  of  the  slope 
jdR/dz'l  «  1.  The  depth  H  is  assumed  finite,  and  it  follows  that  R  is 
small  of  the  same  order,  us  dR/dz'.  In  the  analysis  to  follow  we  shall 
also  require  that  R  be  small  compared  to  the  wavelength  of  the  incident 
wave  system,  or  that  KR  «  1. 

We  wish  to  obtain  the  velocity  potential  of  leading  order  in  the  small 
parameters  of  slenderness  and  oscillation  amplitudes  in  order  to  obtain  a 
consistent  set  of  linearized  equations  of  motion  for  the  body.  However, 
it  will  turn  out  that  the  potentials  of  different  phases  of  the  motion  are  of 
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different  orders  of  magnitude  with  respect  to  the  slenderness  parameter. 
For  example,  the  potential  due  to  surge  or  pitch  is  of  order  R  as  R  — ^  0 
whereas  the  potential  due  to  heave  is  O(R^).  Similar  differences  will  oc 
cur  in  considering  the  components  of  each  potential  which  arc  in  phase 
and  out  of  phase  with  the  respective  velocities  of  the  body.  In  order  to 
circumvent  these  difficulties  without  unnecessary  higher  order  perturba¬ 
tion  analysis,  we  decompose  the  velocity  potential  in  the  following  form; 


$(x,y.z;t)  =  4)^(x,y,z;t)4  4>^(x,y,z;  t)+  (j)^(x,y,z;t) 
•f  A  1  g/o  cos  (Kx  -  u)t)  +  <|>^  (x,  y,  z;  t)] 

O 


[6] 


where  (j)^^  a^re  linear  in  the  displacements  ( ^ ^  ,  4^ )  and  their 

time  derivatives,  respectively.  The  potential  A  g/w  cos  (Kx  -  cot) 

represents  the  incident  wave  system  and  the  potential  A  ( x,  y ,  z;  t)  re¬ 
presents  the  diffracted  wave  potential,  corresponding  to  waves  incident 
on  a  restrained  body.  Each  potential  4^  in  Equation  |  6]  must  sati.'jfy  the 
free  surface  boundary  condition  and  the  radiation  condition;  the  complete 
potential  $  mu^^l  satisfy  the  boundary  condition  on  the  body,  This  condi¬ 
tion,  Equation  |  3]  ,  is  r^educed  as  follows; 


8r'  i)y'  ^ 

dy'  at  f)z'  at 


dR  dz' 
d z ’  at 


a<i>  a^>  dR 

Dr'  a  z '  d  z ’ 


on  r'  -  R(z'), 


or  neglecting  second-order  terms  in  A,  ^^>3  4'> 


a$ 


—  —  -  i  ^  +  (z  -  ^  +  (t  -  xcV)dR/dz  =  0 


on  r  =  R  (  z ), 


|7] 


where  a  dot  denotes  differentiation  with  respect  to  time.  Substituting 
Equation  [6]  into  Equation  [7]  and  separating  terms  according  to  their 
dependence  on  different  displacements,  we  obtain  the  following  boundary 
conditions  on  the  body: 


.  /  94*0  \ 

—  =  I  COS0  +  OiR  - j 

r  V  3z  / 


’‘t’vh  .  /  "'PdjV 

._  =  4,(,_,G)cose  +  o(R— ) 

r0R  /  9<t>r, 

—  =  -  - +  0(R  - ^  ) 

3r  8z  \  0z  / 

^  r  / 

j -  =  —  oj  1  cos  6  sin  wt  —  (  KR 


+  OIR'^) 


9r  8z 


=  -  OJ  fcosS  sinojt  —  f  KR  c  os^  0  H - 


j  c  os  ojtj 


+  O(R^) 


-  OJ  cos  0  sin  ojt  +  ^ ;|  KR  +  |  KR  cos  Zoj  cos  ojtj 

/  ^‘t'A  \  2 

+  0(r— ).0(r2) 

o 

O 

To  satisfy  the  above  boundary  conditions,  we  employ  slender-body 
theory.^  For  example,  the  potential  satisfying  Equation  [8]  is  an  axial 
line  of  horizontal  dipoles,  of  moment  density  ^  ^  [R(z)]^  per  unit  length. 
Thus  in  an  infinite  fluid, 

<!>{:  =  I  4  J  |R(zi)]^  +  (z  -  ^  dzj  [12] 

To  satisfy  the  free  surface  and  radiation  conditions,  we  substitute 

o  y  —  ^ 

for  the  source  potential  [r'^  +  (z  -  Zj)"^]  the  potential  of  an  oscillating 
source  under  a  free  surface.^  With  this  substitution  wc  obtain,,  in  place 
of  Equation  [12]: 

=  H  J  ^  1R(z]  )]^ -^  >^1  r^  +  (z  -  zj  )^]  ^ 

^  i”  ^  Jp  (hr)  dk^  dz^  |l3]i 

o  ®  ° 

+  TTojKfJ^  [RCzi)]^  dz;^ 
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and,  in  a  similar  fashion. 


=  i  +  Jr  +  (2;  -  Zj)^] 

+  TTLlKlj)  J  [R  (Zj^)]^  (z^  -  Zq) 


♦  ^  =  U  J^RU,)f^{|r^+(--..l)h'" 

+  ttcoK^  f  R(z,)  ^1^  J  (Kr)dz 

J.H  dzi  «  ^ 

+  R^  sincjt  “  KR'^  cos  wt  — — V-/ 1  r^  +  ( z 

9^  9x2  J  V 

X  r^k  +  K  k(z+Z|)  ,,  .  ,,  1  , 

'^J„irrr'’  i' j„(kr)dk|<j»i 

-ffcoK  f  KR  +  ~)  R  sinul 

■’J-H  \'  tiz^y 

-  r2  cos  ujt  +  i  KR”^  sinujt  R  (Kr )  dzi 

9x2/  ° 


KR  +  — )  R  coswt 


+  r2  sincjt  “  KR'^  cos  wt 

9x 


r‘^  +  (z  -  z,)  ] 


-  TT  OJ  K 


-  r2  cos  ujt  i  KR”^  sinujt  — — 

9x2 


(Kr )  dzi 


where  ^  denotes  the  Cauchy  principal  value.  From  the  Appendix  we  see 
that  the  potentials  [13]  to  [16]  satisfy  the  boundary  conditions  [8]  to  [11], 
respectively,  with  a  maximum  fractional  error  of  order  R.  Unfortunate¬ 
ly,  this  error  is  not  so  small  as  in  the  classical  slender-body  theory  for 
an  infinite  fluid,  where  the  error  is  of  order  R  log  R;  for  this  reason 

0 

the  present  theory  may  not  hold  for  as  wide  a  range  of  slenderness  as  in 
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the  aerodynamic  case.  However,  for  the  slender  floating  bodies  which 
are  envisaged  at  present  (viz.,  a  rocket  vehicle  or  one  support  of  a  stable 
platform),  this  is  not  expected  to  cause  practical  problems. 

The  values  of  the  potentials  [13]  to  [  16]  on  the  body  may  be  found 
by  setting  r  =  R(z)  and  retaining  the  leading  terms  for  small  R.  To  lead- 
ing  order,  only  the  singular  term  [r  +  (z  -  ^  contributes  to  the  in¬ 

tegrals  over  Z]^,  and  the  integrals  may  be  evaluated  directly  since  for 
any  continuous  bounded  funcliuu  f(zj)  and  small  values  of  r, 

f(Zj)[r'^  +  (z  -  Zj)  ]  ^  dzj  =  -2f(z)logr  +  0(1) 

f  ( Zj)  ~  [r^  +  (z  ~  ^\)^]  ~  cos  0  +  0(1) 

f(z,)  +  (z  -  z,)^]  ^  dzi  =  2  cos  20  +  0(1) 

8x2 

for  -  H  <  z  <  0,  r  «  H. 


I 

I 

I' 


Thus  on  the  body, 

=:  1  R(z)  COS0  +  O(R^) 

[17] 

(+).  "  —  ifi  4-  0(14^^ 

'i|j  'r  \  \  ~  ^ 

[18] 

4)^  -  -i  R  log  R  +  0(r2) 

[19] 

=  CO  e^^  [  ( 1  KR  +  )  R  log  R  cos  cot  +  R  cos  O 

sin  cot]  +  0  (r2) 

=  CO  R  COS0  sin  cot  +  0(r2  log  R) 

[20] 

8 


THE  FIRST-ORDER  FORCES  AND  EQUATIONS  OF  MOTION 


From  Bernoulli's  equation,  the  linearized  pressure  on  the  body  is 


9$ 

P  =  “  P  ~  Pg2 


_  pgz  -  p 


9<t>j^  a<t)j^  P  a(j5^ 


dt 


'  I  T  r  p  *1 

~  P  ~  sin(KR  cos  0  -  tot)J 


pgz  +  p4R(z)cos0  +  pi)jR(z)(z  -  Zq)cos0  -  pAto^e^^R  cosO  cos  cot 
+  pgAe^^oSin  cot  -  pgAKe^^'R  cos  0  cos  cot  +  0(R^  log  R) 

O 

••  ..  lKf7 

Pg2  +  (z)cos0  +  p4'R(z)(z  “  Zq)cos0  +  pgAe  sin  cot 

-  2pco2ao^^R  cosO  cos  ojt  +  0(r2  log  R)  [21] 


The  force  and  moment  exerted  on  the  body  by  the  fluid  are  obtained 
by  integrating  the  pressure  over  the  surface.  In  the  absence  of  any  other 
external  forces,  the  force  or  moment  must  equal  the  respective  accelera.^ 
tion  times  the  mass  or  moment  of  inertia  of  the  body.  Thus,  with  ^  the 
unit  normal  vector  into  the  body,  the  equations  of  motion  are 

m^  =  ff  pcos(n,x)dS  ["’2] 

m(^  +  g)  =  //  p  cos  (n,  z)dS  |  23] 

I  =  //  P  1  ~  zq)  cos  (n,  x»)  -  X  cos  (n,  z)]  dS  [24] 


where  m  is  the  body's  mass,  I  its  moment  of  inertia  about  the  center  of 
gravity,  and  the  surface  integrals  are  over  the  submerged  surface  of  the 
body. 

In  computing  the  pressure  integrals  over  the  body  surface,  it  is  ex¬ 
pedient  to  employ  the  (x',y',z')  system,  fixed  in  the  body.  The  direction 
cosines  are 
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cos(n,x')  =  -  cos  6  +  O(R^) 


cos(n,z')  =  +  0{R^) 

dz' 


and  the  forces  along  the  (x,z)  axis  are  related  to  the  forces  along  the 
(x' ,  z' )  axis  by 


F  X  ^  I  cos  i|)  +  F  |v>sin  i|i  =  F  ,  +  F  ,  +  O(Jj^) 

X  X  Ti  X  '  Z 


F^.  sin  \\f  -  F  -  \\t  F^\  +  0  (i|j  ) 
Thus  the^  equations  of  motion  may  be  written  in  the  form 


me  =  J 

f"j 

'o 

r^*-;  +  x'4- 

cos  6  4 

fr)  pRd.'dO’ 

C 

m(C  +  g)  J 

fZlT 

L  J 

.^*-^.4  x'^ 

-H 

/dR  ^  , 

^^44^cos 

0  j  p.R  dz'  do' 

U\> 


Zir  ;  t  x> 


Jq  J- 

--n 

''0  •'( 


0  J-H 

Zn  r  ^  +  x'i|j 


[{z'  -  Zp)  cos(n,x')  -  x'  cos(n,z')]pRdz'd0 


0 


(z'-  zq)cos  0  pRdz'dO'  4  0(R'^) 


where  is  the  free  surface  elevation  at  the  body.  Substituting  Equation 
(21]  for  the  pressure  and  neglecting  second-order  terms  in  the  oscilla¬ 
tory  displacements  vjj  ,  and  A,  we  obtain 


^  2  TT  ^0  i 

Jn  J  tr  ' 


m^ 


Pg 


cos  0'  4  ijj 


dR 


'0  •'-H 

Zt!  ,.0 


dz 


r)(^' 


4  ;  -  i|jR  COS  0')R  dz'  d0' 


^  CTT  ^  I 

-  p  I  I  COS  0'  R  COS  0'  4  l|jR(z  -  Zq)  COS  0' 
•'0  •'-H 


K  T '  7  K7 ’ 

4  g  Ae^  sin  u)t  —  2co  Ae  R  cos  0'  cos  wt]  R  dz'  d0' 


-  JTpg  j  +  Zvpz' 


p  IT 


r  [|  +  4i(z' -  Zq)  -  2  cos  (.Jt]  R 

J-H 


or,  since 


r„( 


4jR  +  241Z'  ~jR  dz'  -4.  J  ^(R2z')dz 


dz 


-H 


dz ' 


it  follows  that 


me 


=  -  p  r  [e  +  ^(2  ~  ^r)  ~  Zco^Ae^^  cos  wt]  S(z)  dz  + 

J-H 


where 


S(z)  =  TT  [R(z)]‘ 


is  the  sectional  area  function. 

In  a  similar  manner  we  obtain 


m  (C  +  g)  =  -  P  g  ^  S(0) +‘’pg  f  S(z)  dz  +  pg  A  sin  uit  f 

J-H 


+  0(R  log  R) 


“0 

141  --Pg4'  (z  -  Z(~)  S(z)dz 

J-H 


f  +4’  (2  “  Zp)  —  2  cj^Ae^^  cos  wt]  (z  -  z 

-H 


dz' 


'  =  0 


0(R^  log  R) 
[25] 


dz 

|26] 


q) S (z)dz 


+  0(R^  log  R) 


[27] 


From  Archimedes'  principle,  or  equivalently,  satisfying  Equation 
[26]  to  zero  order  in 

gm  =  pg  S(z)dz  [28] 

J-H 

and  thus 

-0 

=  "pguS(0)+pgAsintot  I  dz  +  0  (R"^  log  R)  129] 

J.H 


while,  from  Equations  [28]  and  [25], 

0 

2m^:=-pr  |[[i(z  -  Zq)  -  2oo^Ae'^^  cos  cot]  S  (  z )  dz  +  0  ( R  ^  log  R)  [30] 


Let  us  denote; 

I  =  mk^ 


y  = 

m 

Vertical  Prismatic 

Coefficient 

A 

pHS(0> 

Pn  = 

m 

f 

•'-H 

(z  - 

Zq)"  S(z)  dz 

(n  =1,2) 

rO 

O.JK)  = 

JP_ 

m 

eK== 

(z  -  Zq)"  S(z)dz 

(n  =  0,  1) 

and  note  that  ^  ,  and  4*  rnust  be  sinusoidal  with  frequency  w.  The 
equations  of  motion  then  become 

[31] 

[32] 

[33] 


(1  -  XKH)C  A(1  -  XKHQq  )  sin  wt 
2  ^1  \ 

p2  +  ky  — — j  4*  +  Pj  ^  -  -  2AQj^  cos  ut 
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Note  also  that  surge  and  pitch  are  coupled,  unless  Pj  =  0  or  unless 
the  centers  of  gravity  and  buoyancy  coincide. 

The  above  equations  of  motion  are  not  unexpected.  The  restoring 
forces  on  the  left-hand  side  consist  of  hydrostatic  and  inertial  forces  plus 
entrained  inass  terms  which  double  the  inertial  force  at  each  section. 

This  might  have  been  deduced  as  a  consequence  of  slender-body  theory 
and  the  fact  that  the  entrained  mass  of  a  circular  cylinder  in  an  infinite 

o  o 

fluid  is  just  equal  to  the  displaced  ma  s.s  .  Tn  other  words,  the  hydrodynam¬ 
ic  forces  on  the  left-hand  side  of  Equations  [31]  to  [33]  could  have  been 
obtained  by  neglecting  the  presence  of  the  free  surface.  Moreov^er,  the 
exciting  forces  on  the  right-hand  side  of  these  equations  are  those  which 
follow  from  the  "Fr oude-Krylov"  hypothesis  that  the  pressure  in  the  wave 
system  is  not  affected  by  the  presence  of  the  body.  These  results  are,  of 
co.urse,  a  consequence  of  the  fact  that  the  body  is  slender. 

The  solutions  of  Equations  [31]  to  [33]  are 


^  =  A  sin  tot 


1  -  XQqKH 

1  -Xkh 


34] 


^  =  2  A  cos  tot 


PlQl-QptPz  +  ^y-Pl/K) 

2(p2  +  '"y  -  Pf/K)  -  Pf 


35] 


t|j  =  2  A  cos  tot 


PiQq- 


L  2(P2  +  kj  -  Pj/K)  -  Pf 


36] 


We  note  that  when 


K  = 


XH 


there  is  resonance  in  heave,  and  when 


K  = 


37] 


[38] 
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there  is  r-';sonance  in  pitch  and  surge.  To  determine  the  oscillat'on  am¬ 
plitudes  in  the  vicinity  of  these  resonance  frequencies,  it  is  necessary  to 
consider  the  damping  mechanism  due  to  energy  dissipation  in  outgoing 
waves.  Thus,  for  these  frequencies,  we  must  consider  the  free-surface 
effects  on  the  restoring  forces.  For  this  purpose  we  must  retain  some 
terms  which  are  of  second  ordei  in  the  radius  of  tlie  body. 

THK  DAMPING  FORCES 

The  damping  forces  will  follow  by  considering  the  last  terms  in 
Equations  [13]  to  [16]  and  will  consequently  be  of  higher  order  in  R  than 
those  terms  which  we  retained  in  the  previous  analysis.  This  procedure 
is  nevertheless  consistent,  since  at  resonance  the  lower  order  restoring 
forces  vanish.  In  other  words,  we  are  retaining  the  lowest  order  force 
or  moment  of  each  phase  separately.  For  a  further  discussion  of  this 
point,  see  Reference  4. 

We  proceed,  therefore,  to  study  the  damping  forces,  or  the  forces 
in  phase  with  each  velocity.  The  only  contribution  from  Equations  [  13] 
to  [16]  is  the  potential 

o 

=  irciKe^^  Jo(Kr)dzi 

[39] 

Since  Jq  ( Kr )  =  1  -  j  (Kr  )^  +  .  .  .  ,  it  follows  that  on  the  surface  r  =  R  ( z), 
the  leading  terms  are 

Jq  (Kr)  ^  1 

a 

and 

(Kr)  ^  -  i  K^x  =  -  4  K^Rocos  e 
dx  0  ^  ^ 

Thus,  to  second  order  in  R  the  damping  potential  on  the  body  is 
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4>*  =  TTU)  K 


Kz 


dR 
dz , 


-  \  [4+4j  (Zj  -  Z(-;)]  k2r(z)R^{Z;^)  0  0  3  0^-  dZj 

1  Kz  »  r  K.Z|  dS  , 

=  icoKe^^^  e  1  - —  dzi 

J-H 


[40] 


1  . .  tr  3 


0  y- 

R(z)  cos  0  (  [^  +  4j(z,  -  Zq)]  e  1  S(zj)  dz 

J-H 


The  damping  pressure  on  the  body  is 


-  P 


34-*  ., 
at 


=  -  i wpKe^^i  r 
J-H 


Kz ,  dS 


1  —  dz, 
dz. 


[41] 


+  ■TwpK^R(z)  cos  0  f  [^  +  4^(zi-z;Q)]e  lS(Zj)dzj 

J-H 


Then  the  heave  damping  force  is 


^z  =  p*R|Jd0dz  .  -i..pKi^J^^ 

Similarly,  the  surge  damping  force  and  the  pitch  damping  moment  are 

O 

,0  ^Zir 


[  J 

J-H  J  0 


cos  0  R  d0  dz 


[43] 


)pK^^J  Stz)e^^dzj^J  [^  +  4<(zj  -  zq)]  e^'^^1  S(zj)dz 


and 
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0  .  Ztt 


-  j  j  ~  z^i  cos  G  R  d0  dz 

’-II  Jq 

-  I  wp  ^  J 


[z  -  Zp)  S(z)  e  dz 


^  J  [4  +  4j(Zi  -  Zq)]  S(z^)  dz^ 


or  in  terms  of  the  integrals  Pj  ,  P2,  Qq,  and  Q 


P 


=  -  |copKi[K  ^  Q(^(K)  -  S(0)] 


ojm^K  . 


1  win  ^  . 

XKH]‘ 


i  ^  K^Qq(K)[4Qq(K)  +  ^Q,(K)] 


[441 


[45] 


[46] 


M*  =  -  i  K^Q^{K)[iCiQ(K)  +  ^Qj{K)] 


[47] 


In  place  of  Equations  [31],  [32],  and  [33],  we  obtain  the  damped  equa¬ 
tions  of  motion 


2^  +  Pjv|j  =  -2AQQ(K)cosojt 


+  I  ^  K^Qq(K)  leQo(K)  +  +Qi{K)] 


[48] 


(1-XKH)4  =  a  sin  cot  [  1  -A'KH  Qq  (K)] 

+  i  —  [1  -  Qo(K)XKH]‘ 

'^P  X*h2  ° 


[49] 


.(P2  -  P^/K  +  ky)  4  P|  &  =  -  2A  Qj,(K)  cos 


tot 


+  i  ^  K^QjCK)  [|Qq(K)  +  4jQi(K)1 


[50] 
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The  damping  terms  of  these  equations  of  motion  are  given  by  the  terms 
linear  in  the  velocities  ^  >  cS-Jid  ijj.  It  should  be  noted  that  for  a  slender 

body  rn  — *  0,  and  thus  the  damping  coefficients  will  be  small,  as  was  to  be 
expected.  To  solve  these  equations  for  the  three  unknown  displacements 
and  their  phases  is  a  straightforward  but  tedious  matter.  For  applications 
in  ranges  not  including  a  resonance  frequency,  it  is  much  simpler  to  em¬ 
ploy  the  undamped  equations  of  motion,  [31]  to  [33]  ,  and  the  resulting 
di  spla c ement.s ,  [34]  to  [36]. 


CALCULATIONS  FOR  THE  CIRCULAR  CYLINDER 

As  a  special  case,  we  shall  consider  the  circular  cylinder  R(z)  =  R 
constant.  Then 


P.  = 


Pi  = 


Qo(K)  = 


1 , 0 
1 

~  (z  -  ZQ)dz  -  -  -  Zq 

^  J-H 

^  7  ? 

t  J  (z  -  z^rdz  =  +  Hz^  + 


-H 

,0 

-H 

.0 

-H 


KH 


1  Kz,  ,,  1  -KH  1  ,,  -KHw,  , 

Ql  H  J  ®  ~  ^  e  -  (1  -  e  )( 1  +  Kz^) 


K^H 


We  shall  assume,  moreover,  that  the  centers  of  buoyancy  and  gravity  co¬ 
incide,  or  Zq  =  -  H/2,  so  that  the  equations  of  motion  are  uncoupled  and 
there  is  no  resonance  in  pitch  or  surge. 

Then 


Qj  (K) 


-KH 


-KH 


2  K 


K^H 
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Figure  2  -  Plot  of  the  Surge  Amplitude -Wave  Amplitude 
Ratio  for  the  Circular  Cylinder  .  ' 


Figure  3  -  Plot  of  the  Pitch  Amplitude- Wave  Slope 
Ratio  for  the  Circular  Cylinder 
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and  it  follows  that 


A 

2KH 


( 1  -  e  )  cos  cot 


[51] 


Plots  of  the  above  amplitudes  and  the  heave  phase  angle  are  shown 
in  Figures  Z  to  6  as  functions  of  KH.  Figure  2  shows  the  ratio  of  surge 
amplitude  to  wave  amplitude.  For  zero  frequency  this  ratio  is  one  and 
for  increasing  frequencies  it  decreases  monotonically  to  zero.  Figure  3 
shows  the  ratio  of  pitch  angle  to  the  maximum  wave  slope  KA,  multiplied 
by  the  coefficient  C  =  \  +  6(ky./H)^,  This  coefficient  is  equal  to  one  if  the 
mass  in  the  cylinder  is  uniformly  distributed  throughout  its  submerged 
length.  The  ratio  starts  at  one  for  zero  frequency  and  decreases  mono¬ 
tonically  to  zero.  Thus  the  pitch  amplitude  is  always  less  than  the  wave 
slope.  Figure  4  shows  the  ratio  of  heave  amplitude  to  wave  height  for 
frequencies  away  from  the  vicinity  of  resonance.  Near  resonance,  the 
amplitude  is  shown  in  Figure  5  and  the  phase  angle  in  Figure  6  for  the 
particular  case  R/H  =  0.  1.  The  ratio  of  heave  amplitude  to  wave  ampli¬ 
tude  is  unity  for  zero  frequency,  rises  to  a  maximum  of 

at  the  resonance  frequency  KH  =  1,  and  then  decreases  monotonically  to 
zero.  The  phase  angle  is  similar  to  conventional  one-degree-of-freedom 
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Figure  6  -  Plot  of  the  Heave  Phase  Lag  for  the  Damped 
Circular  Cylinder  with  R/H  =0.1 


harmonic  oscillators  with  linear  damping;  for  low  frequencies  the  heave 
displacement  and  wave  height  are  in  phase,  at  resonance  they  are  in  quad¬ 
rature,  and  at  high  frequencies  they  are  180  deg  out  of  phase. 


DISCUSSION  AND  CONCLUSIONS 

The  damped  equations  of  motion  as  given  by  Equations  [48]  to  [50] 
may  be  solved  for  an  arbitrary  body  of  revolution  to  obtain  the  oscillation 
amplitudes  and  phases.  Except  in  the  vicinity  of  the  resonance  frequencies 
defined  by  Equations  [37]  and  [38],  it  should  be  sufficient  to  use  the  sim¬ 
pler  undamped  equations;  the  resulting  oscillations  are  given  by  Equations 
[34]  to  [36].  Plots  of  these  oscillations  are  shown  in  Figures  2  to  6  for 
a  circular  cylinder,  with  the  important  restriction  that  the  centers  of  buoy¬ 
ancy  and  gravity  coincide.  If  this  restriction  is  relaxed,  a  resonance  will 
be  introduced  into  the  equations  for  pitch  and  surge,  but  the  frequency  of 
this  resonance  may  be  kept  small  by  ballasting.  The  amplitudes  at  reso¬ 
nance  are  extreme,  but  the  resonance  frequency  for  heave  is  quite  small 
and  can  be  kept  out  of  the  practical  range  of  ocean  waves  by  making  the 


21 


draft  sufficiently  large.  It  would  seem  wise  to  do  this  in  practice  and  to 
provide  appropriate  ballast  so  that  the  pitch  resonance  occxirs  at  or  below 
the  heave  resonance  frequency.  From  Equations  [37]  and  [38]  this  re¬ 
quires  that 


+  kj  _  i  p2  xH 

The  advantage  of  spar-buoy-type  bodies  lies  in  their  very  small 
motions  in  the  higher  frequency  range.  By  proper  design  this  advantage 
may  be  utilized;  thus  very  calm  motions  can  be  expected  in  waves. 
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APPENDIX 


Here  the  potentials  4)^,  <j>^  ,  4>^  ,  and  defined  by  Equations  [i3] 

to  [I6],  are  shown  to  satisfy  the  boundary  conditions  [8]  to  [11],  respec¬ 
tively,  to  leading  order  in  R,  For  this  purpose,  let  us  consider  the  po¬ 
tential 


4^ 


9f(zi,t) 


at 


•^[r^  +  (z  -  Zj  )^] 


+ 


+ 


^k(z+zp  dk^  dzj 

f  (z  ,  t)  ( 2  +  Z])  j 


1 


[54] 


whe;r-e  f(z^,t)  has  sinusoidal  time  dependence  with  circular  frequency  oi. 
By  apipropriate  choice  of  the  function  f,  the  potentials  4)^,  4>^  >  and 
4>^  can  all  be  obtained  from  4)  and  94j/9x.  Thus  it  is  sufficient  to  estab¬ 
lish  that  the  following  conditions  are  satisfied  on  the  body  surface  r  =  R: 


94^  ~ 

9r 


R 


[55] 


9^  ip  ^  cos  0  0 

ar  ax  j^2  at 


[56] 


Employing  an  alternative  form  of  the  source  potential,^  we  write  ip 
in  the  form 


-1-  ^^[r2+  (z-  zj)^]'^  +  [r^-  +  {z+z^)"] 

fCO 

— ^  jQ(kr)  dk^  dzj 


1 

2  T  2 


[57] 
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K  r  f(s.  .t)  jQ{Kr)dz, 

J-H 


(57] 

continued 


=  +  4^2 


where. 


4i 

'^z 


X 

2 


K 


+  (Zj  -  z)^']  ^  +  [r^  +  (z^  +  z)^]  ^  j- 

jQ(kr)  dkdzj 
f(zi,t)  jQ(Kr)dZj 


dzj 


The  potential  corresponds  to  an  axial  distribution  of  simple  sources 
together  with  an  image  distribution  above  the  free  surface  z  =  0.  To 
emphasize  this  fact  we  write  i|J]  the  form 


4]  ~  2 


r” 

3  1  1 

I-H 

,  t)  [  r  +  {z  -  z,  r].  dz 


158] 


From  the  conventional  slender-body  theory  of  aerodynamics,  we  may  ex¬ 
pect  this  potential  to  satisfy  the  boundary  conditions  [  55]  and  [56]  on  the 
body  to  leading  order  in  R.  In  fact,  differentiating  with  aspect  to  r  and 

y 

neglecting  terms  which  are  of  order  R  or  R  cos  0  in  the  neighborhood 
of  the  body  r  =  R,  we  have 


9  i(j 
~ar 


1  r  H  8  7  • 

-  =  -  i  — f(-lzi  .t)r[r‘-  +  (z-zj)  ]  ^  dzj 

J  -J-J  ^  * 


9  , 

I,, 


r  r  +  ( z  -  z  1 


3 

\2  1  2 


dz , 


>91 
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z 


n 


i  ^ 
^  9t 


+  (z  -  ) 


-J-H 


i  ^ 

r  at 


and  similarly 

^^^1  _  cos  B 
9r  ax  1-2  at 


159] 

continii  od 


[60] 


Thus  on  the  body  the  potential  i(j  satisfies  the  conditions  [55]  and  [56]  to 
loading  order  in  R  .  To  establish  that  the  same  is  true  of  i|j,  we  now  show 
that  the  contributions  from  i|)2  and  a4j2/0x  ^re  of  higher  order  in  R. 

Since 


—  jQ{kr)  =  -kj^(kr) 


it  follows  that 


9+2 

9r 


-K  r  —X  ■^^l^  J,  (kr)dkdz, 

J-H  9t  Jo  -  K  ^  ^ 

+  ttcoK^  r  f  (z,  ,  t)  Ji  (Kr)dz, 


[61] 


We  wish  to  show  that 


9+2 

0r 


0(f) 


and 


dr  ax 
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as  R  0,  and  thus  that 
94j2 

9r  9r 

for  R/H  «  1.  From  the  series  expansion  of  the  Bessel  function, 

Jj^(kr)  =  ikr  +  O(k^r^) 

and  thus,  where,  this  expansion  is  permissible  in  Equation  [61],  the  re¬ 
sulting  terms  are  clearly  of  order  fK.  However,  in  the  neighborhood  of 
z  =  0,  the  power  series  expansion  is  not  permissible  in  the  integral  over 
k.  It  follows  that,  in  the  neighborhood  of  r  =  R, 


and 


9^2 

3r  3x 


« 


9^  ijjj 
9r  9x 


9^2 

9r 


K 


9f{0,t) 

9t 


e^^l  J^(kr)dkdzj 


+  O(fR) 


S  -  K 

^  -  K 

=  -  K 


9f(0,t) 

at 

9f(0,t) 

9t 

9f(0,t) 

9t 


Ji  (kr) 

— - dk 

k 


-  0(f) 


Similarly, 


^2 
9r  3x 


Thus,  on\the  body, 


and 


'  ^^^1  \ 

3i"  9x  i  0r  9x  j 

Therefore,  the  potential  4’  satisfies  the  conditions  [55]  and  [56]  with  a 
fractional  error  of  order  R. 
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